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1.  INTRODUCTION  AND  PURPOSE 


The  principal  goal  of  the  research  effort  reported  here  Is 
to  derive  certain  stochastic  attrition  models  from  carefully 
stated  probabilistic  assumptions  on  behavior  of  Individual  com- 
batants and  on  interactions  among  combatants.  Although  several 
processes  derived  appear  to  be  stochastic  analogues  of  the  so- 
called  "Lanchester  equations  of  combat,”  In  this  paper  we  do  not 
deal  with  the  square  law/linear  law  dichotomy  that  has  so  long 
engendered  both  debate  and  confusion.  Our  explicit  identifica- 
tion of  the  assumptions  underlying  a number  of  stochastic  models 
of  combat  attrition  permits  determination  of  the  applicability 
of  a given  model  in  terms  of  those  assumptions.  If  a combat 
situation  is  felt  to  be  well  described  by  the  assumptions,  then 
the  model  is  appropriate.  Hence,  there  Is  no  need  to  maintain 
the  poorly  understood  and  misleading  terminology  and  Interpre- 
tation of  "square  law  processes"  and  "linear  law  processes." 
However,  the  processes  derived  here  we  call  LANCHESTER  pro- 
cesses In  honor  of  Lanchester's  fundamental  contributions  to 
the  subject. 

The  paper  Is  organized  In  the  following  manner.  Section  2 
describes  in  detail  the  assumptions  underlying  a very  general 
stochastic  attrition  process,  and  then  characterizes  that 
process  as  a Markov  process.  In  Section  3 we  obtain  a number 
of  specific  (and.  In  some  cases  previously  known)  stochastic 
attrition  models  as  special  cases  of  the  general  model  of 
Section  2.  Some  of  these  models  are  analogous  to  the  deter- 
ministic Lanchester  models;  such  analogies  are  Indicated  In 
order  that  the  reader  may  relate  our  work  to  previous  work 
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on  mathematical  modeling  of  combat  processes.  Finally, 

Section  4 deals  with  some  possible  extensions  of  the  processes 
described  here,  computational  aspects,  and  some  additional 
problems  of  Interest. 


2. 


THE  MATHEMATICAL  MODEL 


I 

t'- 
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In  this  section  we  derive  from  a set  of  carefully  stated 
and  detailed  probabilistic  assumptions  a very  general  stochas- 
tic model  of  combat  between  heterogeneous  forces.  The  model  is 
based  on  a classification  of  shooting  weapons  into  four  classes; 
the  four  classes  are  the  result  of  two  independent  dichotomies 
of  shooting  weapons.  First,  shooting  weapons  are  classified  in 
terms  of  the  qualitative  nature  of  the  rate  at  which  opposition 
weapons  are  engaged.  Here  "rate"  is  in  the  sense  of  the  infini- 
tesimal generator  of  a continuous  time  Markov  process,  as  de- 
scribed in  more  detail  below.  A particular  type  of  shooting 
weapon  will  be  said  to  have  independent  engagement  initiation 
if  the  mean  total  rate  at  which  it  engages  (fires  upon)  opposi- 
tion weapons  is  independent  of  both  the  numerical  size  and  the 
precise  structure  of  the  opposition  force.  In  other  words,  the 
mean  engagement  rate  of  a weapon  with  independent  engagement 
initiation  depends  only  on  the  type  of  weapon.  The  possibility 
that  the  rate  at  which  opposition  weapons  of  a particular  type 
are  engaged  may  depend  on  the  size  and  structure  of  the  opposi- 
tion force  is  not  excluded;  the  definition  is  based  only  on  the 
total  rate  of  engagement.  On  the  other  hand,  a shooting  weapon 
will  be  said  to  have  proportional  engagement  initiation  if  the 
mean  total  rate  at  which  it  engages  opposition  weapons  in  a 

n 

target  force  y * Is  of  the  form  I d.y  where  y.  is 

X n ^ ^ ^ 

the  number  of  opposition  weapons  of  type  J present,. and  the  dj  are 
nonnegative  constants.  In  particular,  the  rate  at  which  opposi- 
tion weapons  of  a given  type  are  engaged  is  directly  proportional 
to  the  number  of  weapons  of  that  type  currently  surviving. 
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We  assume  that  every  type  of  weapon  falls  into  one  of  the 
two  classes  defined  above.  For  some  types  of  weapons  there  is 
little  difficulty  with  this  assumption.  For  example,  an  artillery 
piece  seems  rather  clearly  to  possess  Independent  engagement 
Initiation,  while  a hand-held  antitank  weapon  possesses  propor- 
tional engagement  initiation.  But  in  other  cases  (tanks,  for 
example)  the  problem  is  more  difficult.  In  particular,  the 
category  to  which  a tank  belongs  might  be  the  result  of  tactical 
decisions  during  the  course  of  a battle,  might  change  over  time, 
and  might  even  depend  on  the  structure  of  the  opposing  force. 

Our  model  is  incapable  of  representing  such  phenomena.  In  gen- 
eral, physical  situations  compatible  with  independent  engagement 
initiation  seem  to  be  of  two  types.  Either  a shooting  weapon 
simply  fires  (at  a probabilistically  constant  rate)  at  an  area 
in  which  opposition  weapons  are  known  or  thought  to  be  located, 
or  a shooting  weapon  is  mobile  and  moves  forward  (or  is  pushed 
backward)  in  such  a manner  as  to  maintain  a constant  (overall) 
rate  of  engagement  with  opposition  forces.  In  this  case  every 
opposition  weapon  encountered  is  engaged.  Weapons  with  propor- 
tional engagement  initiation  are  those  in  which  detection  of  an 
opposition  weapon  must  precede  its  being  engaged  and  in  which 
the  time  required  to  detect  a given  opposition  weapon  is  expo- 
nentially distributed.  For  such  weapons  every  detected  opposi- 
tion weapon  is  engaged.  There  exists,  of  course,  the  possibility 
that  some  weapon  fits  well  into  neither  class;  the  model  is  un- 
able to  accommodate  this. 

The  second  classification  of  shooting  weapons  is  based  on 
the  maximum  number  of  opposition  weapons  that  can  be  destroyed 
in  a single  engagement.  If  that  number  is  one,  the  shooting 
weapon  is  said  to  be  single-kill;  if  the  maximum  is  more  than 
one,  the  shooting  "Wedpon  is  termed  a multiple ^cill  weapon.  As 
before,  this  classification  is  independent  of  the  type  of  target, 
so  the  assumption  is  restrictive  to  the  extent  that  a particular 


type  of  weapon  may  in  reality  be  single-kill  against  some  types 
of  targets  but  multiple-kill  against  others.  In  general,  the 
property  of  being  a single-kill  weapon  appears  to  be  compatible 
with  both  Independent  and  proportional  engagement  initiation, 
while  the  property  of  being  a multiple-kill  weapon  is  somewhat 
more  compatible  with  the  former.  However,  multiple-kill  and 
proportional  engagement  initiation  are  not  incompatible;  con- 
sider a weapon  that  makes  no  shots  until  an  opposition  weapon 
is  detected,  but  then  makes  a shot  that  is  sufficiently  power- 
ful to  destroy  many  opposition  weapons. 

According  to  our  classification  scheme  there  are  four 
qualitative  classes  of  weapons: 

IS  - Weapons  with  independent  engagement  initiation  and 
single  kills  per  engagement; 

PS  - Weapons  with  proportional  engagement  initiation  and 
single  kills  per  engagement; 

IM  - Weapons  with  independent  engagement  initiation  and 
multiple  kills  per  engagement; 

PM  - Weapons  with  proportional  engagement  initiation  and 
multiple  kills  per  engagement. 

We  next  present  the  assumptions  and  notations  for  our 
generalized  attrition  process.  Further  discussion  of  the  basic 
classification  appears  in  the  comments  on  the  assumptions.  The 
combat  involves  two  heterogeneous  forces  (Red  and  Blue);  weapon 
types  are  distinguished  only  into  the  four  qualitative  classes 
described  above  and,  within  each  class,  by  variations  of  numeri- 
cal parameters. 

Assumptions 

1.  The  Blue  force  consists  of  weapon  types  of  the  class  IS, 
weapon  types  of  class  PS,  weapon  types  of  the  class 
IM,  and  weapon  types  of  the  class  PM.  Hence,  the  Blue 
side  has  altogether  M » M^^  + M2  + M^  + M^^  types  of  weapons. 
The  analogous  numbers  for  the  Red  side  are  N^,  N^,  N^, 
and  N » + N2  + Ng  + N^. 
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The  notation  used  below  Is  of  the  fol.  •)wlng  form: 
the  letter  ”1"  is  used  as  a generic  index  for  Blue  weapon 
types.  Types  1,  are  the  IS  weapons;  + 1, 

Ml  + M^  are  PS  weapons;  weapon  types  + M2  + 1,  . . . , M^ 

+ M2  + M^  are  of  class  IM;  and  the  remaining  types  belong 
to  the  class  PM.  Certain  parameters  below  are  defined  only 
for  Blue  weapons  belonging  to  one  class,  e.g.,  if  that  class 
is  IM,  the  relevant  index  1 ranges  over  1,  ...,  M^;  when  all 
Blue  weapon  types  are  considered,  the  index  1 ranges  over 
1,  ...,  M.  For  Red,  a similar  situation  holds  for  the 
index  J . 

Blue  forces  are  denoted  by  vectors  x € N^,  where  N = 

{0,1,2 ,...};  x^  is  the  number  of  type-1  weapons  (1=1,...,M) 

currently  present.  Red  forces  are,  analogously,  denoted  by 
N 

vectors  y € N . 

2.  Times  between  engagements  initiated  by  a surviving  Blue  IS 
type-1  weapon  are  Independent  and  identically  exponentially 
distributed  with  mean  l/r^Cl),  1 » 1,  ...,  M^^. 

3.  When  a Blue  IS  weapon  of  type  i initiates  an  engagement, 
it  attacks  exactly  one  Red  weapon.  If  the  current  Red 
force  is  the  vector  y the  probability  that  the  weapon  at- 
tacked is  of  type  J is  ag(l,y;J);  the  particular  type-J 

rf  weapon  attacked  is  chosen  uniformly.  Both  choices  are 

I independent  of  the  past  history  of  the  process  and  of 

f each  other. 


4.  The  conditional  probability  that  a Blue  IS  type-1  weapon 
kills  a Red  type-J  weapon,  given  an  attack  on  that  weapon, 
is  Pg(i>j)  for  i “ 1,  ...,  M^  and  J * 1,  ...,  N. 

5.  Assumptions  2,  3*  and  H hold  for  Red  IS  weapons  with  param^ 

eters  rg(J),  J ■ 1,  •••»  J " 1*  •••» 

1 ■ 1,  ...,  M and  X € N^,  and  Pj^(J,l),  J ■ 1,  ...»  N^; 
i - 1,  ...,  M. 
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6.  The  time  required  for  a particular  surviving  Blue  PS 
type-1  weapon  to  detect  a particular  surviving  Red  type-J 
weapon  is  exponentially  distributed  with  mean  l/dn(l,J) 

O 

for  i * 1,  ..•»  M2  and  j » 1,  N.  A Blue  PS  weapon 

detects  different  Red  weapons  Independently  of  one  an- 
other and  can  engage  an  opposition  weapon  only  after 
detecting  it. 

7.  A Blue  PS  weapon  engages  every  Red  weapon  it  detects.  The 
conditional  probability  that  a Blue  PS  type-i  weapon  kills 
a Red  type-J  weapon,  given  detection  and  attack,  is  kg(i,J) 

for  1 » 1,  ...,  M2;  J = 1,  ...,  N. 

8.  Red  PS  weapons  satisfy  Assumptions  6 and  7 with  mean  de- 
tection times  l/dj^(J,l)  and  kill  probabilities  kj^(J,l), 
defined  for  J * 1,  ...,  N2  and  1*1,  ...,  M. 

9.  Times  between  engagements  initiated  by  a surviving  Blue 
IM  type-1  weapon  (1=1,...,M^)  are  independent  and  identi- 
cally exponentially  distributed  with  mean  l/r*(l). 

10.  Given  that  a Blue  IM  type-1  weapon  initiates  an  engage- 
ment against  a currently  surviving  Red  force  y,  the 
probability  that  the  surviving  target  force  has  composi- 
tion z is  ^g(l,y;z).  The  numbers  of  surviving  Red  weapons 
are  otherwise  independent  of  the  past  history  of  the 
attrition  process.  Symbolically, 

^o(l»yjz)  * P{ surviving  Red  force  is  z| engagement 
“ initiated  by  Blue  IM  type-i  weapon 

against  Red  force  with  composition  y}. 

Here,  y,z  € N^. 

11.  Red  IM  weapons  obey  Assumptions  9 and  10,  with  mean  engage- 
ment times  l/rt(J),  J ■ 1,  ...,  N_,  and  kill  distributions 

A M 

^j^(J»3t;w)  defined  for  J ■ 1,  ...,  N^;  x,w  € N . 
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12.  The  time  required  for  a particular  surviving  Blue  PM 
type-i  weapon  (1=1, . . . ,M^)  to  detect  a particular  sur- 
viving Red  type-J  weapon  (J»1,...,N)  is  exponentially 
distributed  with  mean  l/d*(l,J).  A Blue  PM  weapon  de- 
tects different  Red  weapons  Independently  of  one  another 

13.  At  the  Instant  of  each  detection,  a Blue  PM  weapon 
initiates  an  engagement. 

14.  Each  such  engagement  is  Independent  of  the  previous 
history  of  the  attrition  process.  If  a Blue  PM  type-1 
weapon  engages  a Red  force  of  composition  y after  having 
detected  a Red  type-J  weapon  (J=1,...,N),  the  probabil- 
ity that  the  Red  force  surviving  the  engagement  is  z is 
denoted  by  yg(l,J,y;z). 

15.  Red  weapons  of  the  class  PM  satisfy  Assumptions  12,  13, 

and  14  with  parameters  dl(J,l)  and  m„(J,1,x;w)  defined 

^ ^ M 

for  J » 1,  ...,  N^;  1 = 1,  ...,  M,  and  x,w  € N . 

16.  The  detection,  engagement,  and  kill  processes  of  all 
weapons  are  mutually  independent . 

To  avoid  unnecessary  proliferation  of  notation,  certain 
parameters  are  denoted  by  an  asterisk  for  multiple-kill  weap- 
ons and  no  asterisk  for  single-kill  weapons. 

Comments  on  the  Assumptions 

It  follows  from  Assumptions  2,  6,  9«  12,  and  I6  that  the 
time  required  for  a given  shooting  weapon  to  initiate  an  en- 
gagement is,  conditioned  on  survival  of  the  weapon  until  that 
time,  exponentially  distributed.  The  parameter  of  that  expo- 
nential distribution  depends  on  the  class  and  type  of  the 
shooting  weapon,  but  is  one  of  two  forms  according  as  the 
weapon  has  Independent  or  proportional  engagement  initiation. 
Table  1 illustrates  for  single-shot  weapons  of  type  1 on  the 
Blue  side  against  a Red  force  y.  The  same  holds  for  multiple 
shot  weapons  but  with  parameters  having  asterisks. 
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Table  1.  EXPONENTIAL  PARAMETERS  FOR 
OVERALL  ENGAGEMENT  RATES 


Form  of  Engagement 
Initiation  

Parameter 

Independent 

rgd) 

Proportional 

1 dg( i , j )yj 

Recalling  that  the  parameter  of  an  exponential  distribu- 
tion is  the  inverse  of  the  expectation  and  that  the  latter  is 
in  units  of  time,  we  may  interpret  the  entries  in  Table  1 as 
expected  overall  rates  of  engagement  of  the  Red  force  by  the 
weapon  in  question.  Expected  rates  of  engagement  of  particu- 
lar types  of  Red  weapons  are  given  in  Table  2. 


Table  2.  ENGAGEMENT 
OPPOSITION 

RATES  FOR  PARTICULAR 
WEAPON  TYPES 

Rate  of  Engagement  of  ■ 

Form  of  Engagement 

Opposition  Weapons  of  1 

Initiation 

^ Type  j 

Independent 

rg(i)  ag(i,y;j)  j 

Proportional 

dg(1.j)yj  ] 

Derivation  of  the  kill  distributions  Pg  and  Pg  is 

a problem  that  we  have  not  yet  considered  in  depth.  In  applica- 
tion of  our  model  to  computerized  simulations,  this  would  be  the 
problem  most  in  need  of  attention.  Here  we  have  indicated  ab- 
stractly functional  dependences  that  seem  plausible  and,  hence « 
those  that  we  feel  can  safely  be  igpiored.  For  example,  the  kill 
distributions  do  not  depend  on  the  structure  of  the  force  to 
which  a shooting  weapon  belongs,  though  in  principle  they  could. 


Indeed,  a reasonable  way  of  describing  such  dependence  might 
allow  representation  of  the  synergistic  Interactions  of  weap- 
ons on  the  same  side.  Similarly,  the  kill  distributions  for 
PM  weapons  indicate  that  the  distribution  of  weapons  killed 
can  depend  on  the  particular  (type  of)  weapon  first  detected; 
differing  ammunition  or  tactics  used  against  different  detected 
weapons  can  thus  be  modeled. 

In  Section  4 we  discuss  at  some  length  certain  plausible 
forms  for  the  fire  allocation  distributions  ag  and  a^.  In  par- 
ticular, all  distributions  satisfying  some  rather  weak  restric- 
tions are  of  a certain  simple  form. 

All  "engagements"  occur  Instantaneously,  with  ensuing 
total  loss  of  contact.  This  is  admittedly  an  unrealistic  fea- 
ture of  the  model,  although  most  other  models  do  not  seem  to 
have  successfully  addressed  the  difficulty  either.  One  way  of 
including  binary  engagements  with  exponentially  distributed 
lengths  is  discussed  in  [8]. 

Results 

Our  notation  and  terminology  concerning  Markov  processes 
are  those  of  [2],  [3]  and  [6];  the  reader  is  referred  to  these 
for  background  material.  We  recall  here  the  following  facts. 
Let  (X^)  be  a continuous  time  Markov  process  with  countable 
state  space  and  no  Instantaneous  states.  The  matrix  A defined 
by 

P{X.  - J |X-  - 1}  - 1(1, J) 

A(1,J)  - 11m  5 ^ , 

t*0 

where  I denotes  the  identity  matrix,  is  called  the  infinitesimal 
generator  of  (X^).  If 

Pt(i.J)  - P{X^  • J 1^0  “ 

% 
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denotes  the  transition  function  of  (X^),  then  A satisfies  the 
backward  equation 

^t  “ 

and  the  forward  equation 

P'  - P^A  . (1) 

Matrices  are  differentiated  componentwise  with  respect  to  t. 
From  (1)  it  follows  that 

I Pt(i.‘c)A(kJ) 
k 

for  small  values  of  h,  so  that  A(k,J)  may  be  Interpreted  as 
the  (infinitesimal)  rate  at  which  the  process  tends  to  move 
from  state  k to  state  J.  The  quantity  X(i)  ■ -A(i,i)  is 
finite,  but  possibly  zero,  for  each  i.  When  the  process  enters 
state  i it  remains  there  for  an  exponentially  distributed  ran- 
dom time  that  is  independent  of  the  entire  past  history  of  the 
process  and  has  expectation  X(l)~^.  If  X(l)  » 0,  then,  once 
entered,  state  i cannot  be  left  and  is  called  absorbing.  If 
X(i)  > 0,  then  state  1 is  called  stable;  at  the  end  of  the 
sojourn  the  next  state  to  be  entered  is  state  J with  probabil- 
ity Q(1,J)  » A{1,J)/X(1).  The  sequence  of  states  entered  is 
thus  a Markov  chain,  called  the  embedded  Markov  chain,  with 
transition  matrix  Q.  We  refer  the  reader  to  the  previously 
mentioned  sources  for  further  details. 

We  can  now  describe  and  characterize  the  stochastic  attri- 
tion process  engendered  by  Assumptions  1 through  16  above.  Let 
E . consist  of  states  denoted  by  a ■ (x,y);  a is  to  be 

thought  of  as  a possible  pair  of  surviving  forces  at  some  in- 
stant— with  X describing  the  Blue  force  and  y the  Red.  The 
sample  space  for  our  attrition  process  is  the  family  n of  func- 
tions from  R.  to  E that  are  right-continuous  and  have  limits 
from  the  left  with  respect  to  the  discrete  topology  on  E.  The 
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coordinate  (vector-valued)  stochastic  process  ((B. ,R. (here 
M N V u 

Bt’.Q  N and  N for  each  t)  has  the  usual  Interpretation 

Is  the  surviving  Blue  force  at  time  t;  R^,  the  Red  force  at 

time  t.  We  further  define 

Et  “ 0 < u < t}  , 

which  Is  the  history  of  the  attrition  process  until  time  t,  and 

g - a{(B^,R^);  u > 0}  , 

which  Is  the  entire  history  of  the  process.  For  each  a € E, 
we  denote  by  P®  the  probability  law  on  (fi,g)  of  the  attrition 
process  governed  by  Assumptions  1 through  16,  subject  to  the 
Initial  condition 


p“{(Bq,Rq)  - a}  - 1. 

Here  Is  our  main  result; 

THEOREM.  Subject  to  Assumptions  1 through  16  listed  above, 
the  stochastic  process 

(B,R)  - -^2) 

Is  a Markov  process  with  state  space  E.  The  Infinitesimal 
generator  A of  the  process  Is  of  the  following  form.  Consider 
two  states  o ■ (x,y)  and  a'  ■ (x',y')  such  that  y'  ■ y,  * 

- 1 for  some  1 and  » x^^  for  all  k W 1 (the  new  state  a' 

Is  reached  from  a by  the  destruction  of  exactly  one  Blue  type-1 
weapon);  then 


N. 


A(a,a')  » I Tp(J )ap(J ,x;l)pp(j ,l)y 


J-1 


'J 


N, 


N- 


+ I rS(J  ) <^p(J  ,x;x')y 


J»1 


Nj^+N2+j 


N, 


M 

* k-1  J-1  ^^N^+N2+N^+J  • 

For  any  other  state  a'  » (x',y')  with  y'  » y,  x'  < x^  for  all 

i,  and  I (x  - x')  > 2,  which  corresponds  to  the  simultaneous 
i ^ ^ 

destruction  of  more  than  one  Blue  weapon  and  can  be  effected 
only  by  Red  weapons  of  classes  IM  and  PM,  we  have 


A(a.a')  - I +N2+J 

J ^ 


N. 


M 4 

+ \ d»(J  ,k)Wj^(J  ,k,x;x  )yN^+N2+N^+J  ’ 


(4) 


Similarly,  for  a'  * (x',y')  such  that  x'  » x,  yj  » yj  - 1 
and  y^  » y^^  for  J » 


M, 


A(a,o')  - I rg(l)ag(l,y;J )pg(l,J )x^ 


1-1 


M, 


+ yj  ^*Mj^+l 


M, 


I 

B 


i 

M 
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while  for  any  other  state  o'  » (x 
for  all  J and  X “ yp  1 2, 


y')  such  that  x 


Moreover,  for  all  a'  a and  not  of  the  forms  above,  we 


have 


and,  finally 


We  omit  the  proof  of  the  Theorem,  which  one  constructs  by 
using  the  results  and  by  the  methods  of  the  appendix  of  KARR  [8] 
The  expressions  for  the  generator  A have  probabilistic  Inter- 
pretations and  are  written  In  what  we  feel  is  the  most  revealing 
form.  Consider,  for  example,  the  term  A(a,a')  of  Equation  (3). 
Here  the  first  summand  (of  four)  on  the  right  is  the  instantane- 
ous rate  at  which  Blue  type-1  weapons  are  being  killed  by  Red 
IS  weapons  when  the  two  forces  have  compositions  x and  y,  re- 
spectively, and  the  second  term  Is  a similar  rate  arising  from 
Red  PS  weapons.  For  kills  caused  by  single-shot  weapons,  "rate 
of  kill"  and  "rate  of  kill  one  at  a time"  are  the  same.  This 
Is  not  so  for  multiple-shot  weapons.  Hence,  the  third  summand 
on  the  right-hand  side  of  Equation  (3)  must  be  Interpreted  as 
the  Instantaneous  rate  at  which  Blue  type-1  weapons  are  being 
killed  exactly  one  at  a time  by  Red  IM  weapons.  The  Interpre- 
tation of  the  fourth  summand  Is  then  analogous;  note  that  the 


kill  of  a type-1  weapon  can  arise.  In  principle,  from  the  detec- 
tion of  any  type  of  weapon.  Hence  A(a,a')  Is — when  the  struc- 
ture of  the  two  forces  Is  o » (x,y) — the  Instantaneous  rate  at 
which  Blue  type-1  weapons  are  being  destroyed  precisely  one  at 
a time.  We  emphasize  that  single  kills.  In  general,  can  arise 
from  multiple-shot  weapons,  as  Indicated  by  the  presence  of 
the  third  and  fourth  summands  In  Equation  (3). 

For  the  term  of  A given  by  Equation  (4),  only  multiple-kill 
weapons  need  be  considered,  and  interpretations  are  similar  to 
those  given  above. 

The  expressions  A(a,a")  given  in  the  Theorem  are  valid 
only  If  a Is  not  an  absorbing  state.  In  this  process  the  ab- 
sorbing states  are  those  of  the  forms  (0,y)  and  (x,0),  which 
correspond  to  annihilation  of  the  Blue  and  Red  sides,  respec- 
tively. For  such  states  a,  A(a,a^)  = 0 for  all  a'. 

We  next  list,  for  purposes  of  reference  and  completeness, 
some  consequences  of  the  Theorem. 

COROLLARIES.  (a)  For  the  Markov  process  (B,R),  all  states 
of  the  form  (x,0)  with  x € or  of  the  form  (0,y)  with  y € 
are  absorbing;  all  other  states  are  stable  and  transient; 

(b)  If  a (0,0),  then 

P**{(B^,R^)  ■ (0,0)  for  some  t}  « 0 ; 

(c)  For  any  a,  with  - probability  1,  each  component  of 
the  sample  paths  t -*•  B^  and  t -►  R^  Is  nonincreasing.  (There  Is 
no  provision  for  reinforcements;  one  Indication  of  how  this 
problem  might  be  handled  Is  given  In  [8].) 

For  simplicity  we  have  not  written  down  explicit  expressions 
for  the  jump  rates  X(a)  or  the  transition  matrix  Q of  the  embed- 
ded Markov  chain.  These  may  be  obtained  in  a straightforward 
manner  from  the  Infinitesimal  generator  using  the  formulae  given 
before  the  Theorem. 
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^ The  great  generality  of  the  Theorem  possibly  obscures  Its 

[ usefulness.  In  the  next  section  we  Illustrate  that  usefulness 

^ by  setting  forth  In  some  detail  a number  of  special  cases. 
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3.  SOME  SPECIFIC  LANCHESTER  PROCESSES 


The  objective  of  this  section  Is  to  Illustrate  the  scope 
and  utility  of  the  attrition  process  formulated  In  the  previous 
section  by  deriving  (as  special  cases)  a number  of  simplified 
stochastic  attrition  processes.  Some  of  these  processes  are 
known,  principally  In  the  context  of  stochastic  analogues  of 
the  deterministic  attrition  models  of  F.  W.  LANCHESTER  and 
many  others;  the  reader  Is  referred  to  [5]  and  [7]  for  surveys 
of  this  work.  We  point  out  such  connections  In  this  section 
not  to  Justify  the  work  reported  here  but  to  Indicate  Its 
relation  to  previous  research  In  the  same  general  area.  In 
particular,  we  believe  that  our  assumptions  and  weapon  classi- 
fications should  be  taken  on  their  own  merits  and  Judged 
neither  as  stochastic  analogues  used  to  Justify  Lanchester's 
equations  nor  as  a stochastic  Interpretation  of  the  classical 
"linear  law/square  law"  dichotomy.  Applicability  of  our  as- 
sumptions to  any  given  physical  situation  Is  Independent  of 
such  questions. 

The  reader  should  not  Infer  from  the  discussion  above 
either  that  there  Is  no  relation  between  stochastic  and  deter- 
ministic models  of  attrition  or  that  the  relation  Is  uninteres- 
ting or  unimportant.  In  fact,  the  relation  Is  both  Interesting 
and  Important;  details  of  It  are  explored  In  [12]. 

For  each  process  below  we  Indicate  how  to  obtain  It  as  a 
special  case  of  the  attrition  process  described  In  Section  2, 
present  explicitly  the  infinitesimal  generator  A,  Jump  function 
X and  transition  matrix  Q of  the  embedded  Markov  chain,  and 
discuss  the  relation  of  the  process  to  previously  published 


work.  We  begin  with  the  simplest  case  of  homogeneous  forces, 
in  which  there  is  but  one  type  of  weapon  on  each  side. 

Homogeneous  Processes 

Throughout  this  discussion  we  assume  that  there  is  only  one 
type  of  weapon  on  each  side.  The  attrition  process  is  denoted  by 
( (B^ ,R^ ) , where  is  the  number  of  Blue  weapons  surviving 
at  time  t and  is  the  number  of  Red  weapons  surviving  at  time 
t,  and  has  as  its  state  space  E = N x N.  Rotational  dependences 
of  parameters  of  Section  2 on  weapon  type  are  suppressed  for 
simplicity. 

1 . Independent  Engagement  Initiation.  Single  Kills.  In 
this  case  times  between  engagements  initiated  by  a Blue  weapon 
are  Independent  and  identically  exponentially  distributed  with 
expectation  1/rg  and  fatal  with  probability  p^;  corresponding 
parameters  for  Red  are  r^  and  Pg,  respectively.  Then,  the 
attrition  process  has  Infinitesimal  generator  A given  by 

A((i,J),  (1,J-D)  * irgPg 

A((1,J),  (1,J))  = - (IrgPg+JrgPg) 

A((1,J),  (1-1, J))  = Jrgpg  , 

Jump  function  \ given  by 

X(1,J)  * “^^R^R  * 

and  transition  matrix  Q of  the  embedded  Markov  chain  given  by 

Q((1,J),  (1,J-D)  » IrgPg  + JrgPg 

J 

9(U.J),  (l-l.J))  - ■ 

These  expressions  hold  only  for  i > 0,  J > 0.  All  other  states 
are  absorbing. 
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This  particular  attrition  process  was  proposed  by  R.  N. 
SNOW  [13]  as  a stochastic  analogue  of  Lanchester's  square  law 
of  attrition. 


2.  Proportional  Engagement  Initiation.  Single  Kills.  All 

v/eapons  of  the  single  type  on  each  side  are  of  class  PM.  The 
expected  times  to  make  a detection  and  kill  probabilities  are 
dg^jkg  and  d~^,kpj  for  Blue  and  Red  weapons,  respectively.  In 
this  case  the  attrition  process  ^ ^ ^t>0  infinitesimal 

generator  A given  by  ~ 


A((i,J),  (1,J-D)  = iJdgkg 

A((i,J),  (i,J))  = - iJ(dgkg+dj^kp) 
A((i,J),  (i-l,J))  = IJdpkp  , 

Jump  function  X given  by 


X(i,J)  = ij  (dgkg+dj^kpj)  , 

and  transition  matrix  Q of  the  embedded  Markow  chain  given  by 


Q((i,J),  (IJ-D) 


Q((i,J),  (i-l,J)) 


^B*^B 


The  Infinitesimal  generator  of  this  process  is  similar  to  the 
differential  model  known  as  Lanchester's  linear  law,  so  one 
could  interpret  within  the  Lanchesterean  context  this  process 
as  a stochastic  linear  law  process.  Prom  a computational 
standpoint  the  embedded  Markov  chain  of  this  process  is  particu- 
larly tractable:  it  is  a spatially  homogeneous  random  walk  on 

the  2-dlmenslonal  lattice  N x N;  considerable  theoretical  and 
computational  knowledge  exists  concerning  such  processes.  For 
example,  the  probability  that  a given  side  wins  according  to 
criterion  that  Involves  only  numbers  of  weapons  and  the 
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associated  terminal  distribution  of  the  process  can  be  computed 
using  only  the  embedded  Markov  chain. 


3.  Asymmetric  Engagement  Initiation,  Single  Kills.  The 
single  Blue  weapon  type  Is  of  class  IS  with  parameters  r_,p_; 

D D 

the  single  Red  weapon  type  Is  of  class  PS  with  parameters 
dj^.kj^.  For  the  resultant  Markov  attrition  process  the  Infini- 
tesimal generator  A Is  given  by 

A((1,J),  (1,J-D)  = IrgPg 

A((1,J),  (1,J))  = - (IrgPg  + Ud^kg) 

A((1,J),  (l-lj))  = IJdgkg  , 

the  Jump  function  \ Is  given  by 

X(1,J)  = i^gPg  + IJdgkg  , 


and  the  transition  matrix  Q of  the  embedded  Markov  chain  Is 
given  by  , 


Q((1,J),  (1,J-D) 


^bPb 


Q((1,J),  (1-1, J)) 


JdRkR 


TgPg  + Jdgkg 


Observe  that  Q((1,J),*)  Is  Independent  of  1,  which  leads  to 
simplifications  of  certain  cpmputatlons.  In  Lanchesterean  terms, 
this  process  Is  a stochastic  analogue  of  the  mixed  law  process 
proposed  by  S.  DEITCHMAN  [4]  as  a model  of  guerilla  ambushes. 


4.  Independent  Engagement  Initiation.  Binomially  Distri- 
buted Multiple  Kills.  Here  the  single  tyye  of  weapon  on  each 
side  belongs  to  class  IM;  the  exponential  firing  rates  are 
rg  and  rg  for  the  Blue  and  Red  sides,  respectively.  The  kill 
distributions  of  Assumptions  10  and  11  are  binomial.  That  Is, 


v>g(j;k)  = P{k  Red  weapons  survive  an  engagement  Initiated 
by  a Blue  weapon  at  the  start  of  which  J Red 
weapons  are  surviving} 


where  0 < < 1.  That  is,  each  currently  surviving  Red  weapon 

is  killed  with  probability  q^.  Independently  of  all  the  others. 
We  denote  by  q^^  the  analogously  defined  kill  probability  for 
weapons  on  the  Red  side.  The  Markov  attrition  process  then  has 
Infinitesimal  generator  A given  by 

A((1,J),  (i,0)  - t = 0,  J-1 

A((1,J),  (i,J))  - - [lr|(l-(l-qg)J)  + Jr»(l-(l-qg)^)] 

A((1,J),  (k,J))  - . k - 0,  ...»  1-1  , 

Jump  function  A given  by 

X(1,J)  - lr»(l-(l-qg)'J)  + Jr»(l-(l-qg)^)  , 
and  transition  matrix  Q of  the  embedded  Markov  chain  given  by 


An  attempt  to  classify  this  process  as  a "Lanchester  square" 
or  a "Lanchester  linear"  process  helps  to  Illustrate  not  only  the 
difficulties  but  also  the  Irrelevance  of  the  classification.  On 
the  one  hand  the  form  of  the  Jump  function  A is  similar  to  coef- 
ficients in  the  Lanchester  square  equations,  but  for  small  values  ^ 

of  q„  and  qo  the  principal  terms  of  the  generator  are  those  cor-  i 

D A 

responding  to  one  Red  casualty:  namely. 


\ 


A((1,J),  (1,J-D)  - 

- lJr|qg(l-q3)J-^ 
and  to  one  Blue  casualty:  namely, 

A((i,J),  (i-l,J))  - ijrgqg(l-qg)^"^  . 

These  latter  terms  clearly  resemble  the  Lanchester  linear  equa- 
tions, so  no  clear  choice  Is  possible.  We  again  emphasize  that 
within  our  formulation  there  Is  no  need  for  such  considerations; 
assumptions  are  always  dealt  with  explicitly  and  on  their  own 
merits. 

Heterogeneous  Processes 

In  the  processes  discussed  below  there  are  more  than  one 
type  of  weapon  on  at  least  one  of  the  two  sides.  For  such 
processes  the  parameters  (fire  allocations  for  IS  weapons) 

(survivor  distributions  for  IM  weapons)  and  Vg»Ug  (sur- 
vivor distributions  for  PM  weapons)  must  be  specified;  this  Is 
a problem  on  which  our  research  has  not  concentrated.  We  have 
Implicitly  assumed  that  certain  dependences  can  be  neglected; 
for  example,  none  of  the  above  distributions  depends  on  the 
current  structure  of  the  force  to  which  the  shooting  weapon 
belongs.  Below  we  Indicate  some  possible  specific  forms  for 
the  distributions  required.  We  hope  others  will  contribute 
to  this  effort. 

1 . Independent  Engagement  Initiation,  Uniform  Fire  Alloca- 
tion, and  Single  Kills.  All  weapons  on  each  side  are  of  class 
IS,  but  each  side  may  have  several  types  of  weapons.  The  engage 
raent  rates  are  rg(l),  Tg(J)  for  appropriate  values  of  1 and  J; 
the  kill  probabilities  are  Pg(l,J),  Pg(J,l).  Fire  allocation  Is 
uniform  In  the  numerical  sense  and  independent  of  the  type  of 
shooting  weapon.  That  Is,  for  all  1,J  and  y, 
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I 


I it 

T 


13(1. y;J)  - yj/Z  y^ 


While  for  all  1,J  and  x 


a^CjfX;!)  - x^/Z  . 

The  assumption  of  uniform  fire  allocation  Is  plausible  If 
different  types  of  weapons  are  not  significantly  different  In 
their  characteristics  as  targets,  but  may  differ  In  their  char- 
acteristics as  shooting  weapons.  The  resultant  attrition 
process  has  Infinitesimal  generator  A given  by 

y,  ”1 

A((x,y);  (x;y,  . ,y.-l ))  - — i-  I dnCDPod.j  )x. 

^ ^ ”1  Z y,  1-1  ^ 


M,  N, 


"1  lty.d„(l)p„(l,J)x. 
U(x,y);  (x.y))  - - Z Z ^ 

1-1  J-l/  y y. 


Xidp(J )pp(j ,l)y. 


X '1 

A((x,y);  (x^ x -l,....,x„  ;y))  » I dp(J  )p„(J  ,l)y  . , 

1 Z 


Jump  function  X given  by 


“1  "1  ( 

X(x,y)  • z Z 

1-1  J-1  ^ 


dg(l)Pg(l,J)  ^ dg(J)Pg(J,l)l 


and  transition  matrix  Q of  the  embedded  Markov  chain  given  by 


$ 


M, 


Q((x,y);  (x;yj^,. . . . . ,y^^)) 


yj  dgCDPgd  J )Xj^ 


A(x,y)  I y- 


Q((Xjy)j  >X|^  jy ) ) 


’‘i  ‘^R<'>>PR<J-^>yj 

^(x,y)  I Xj^ 


We  do  not  attempt  a Lanchesterean  classification  of  this  process. 

Some  other  methods  of  fire  allocation  are  discussed  In 
Section  4 below. 


2.  Proportional  Engagement  Initiation,  Single  Kills. 

Weapons  on  each  side  are  all  of  class  PS,  with  expected  times  to 
detection  denoted  by  d“^(l,J),  d“^(J,l)  and  kill  probabilities 
denoted  by  kg(l,J),  kj^(J,l).  The  ranges  of  the  Indices  are 
1*1,  . . . , Mj  and  J ■ 1,  ...,  N^.  The  attrition  process  that 
results  Is  a straightforward  extension  of  Process  2 In  the  sub- 
section on  homogeneous  processes.  The  Infinitesimal  generator  A 
Is  given  by 


M. 


^((x,y);  (x;y^^,. . . ,yj-l,. . . ,yj^^))  - y^  dg(l,J  )kg(l,J 


”2  ^2 


i((x,y);  (x,y))  • - I I X y.[d„(l,J)k„(l,J) 

1-1  J-1  ^ J ® ® 


+ dj^(J,l)kj^(J,l)] 


N, 


A((x,y);  (Xj^,...,Xj^-l,...,x  ;y))  • I f »i)yj  » 

t J — 1 

and  other  parameters  may  be  computed  In  the  manner  Indicated  in 
Section  2;^  For  the  sake  of  brevity  we  omit  further  details. 
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3.  Independent  Engagement  Initiation,  Uniform  Fire  Alloca- 
tion, Binomially  Distributed  Kills  of  a Single  Weapon  Type.  This 
example  Illustrates  how,  for  weapons  of  class  IM  In  the  hetero- 
geneous case,  the  kill  distributions  can  be  constructed  so 

D K 

as  to  represent  the  effects  of  both  allocation  of  fire  and  kill 
probabilities.  We  assume  here  that  all  weapons  are  of  class  IM, 
that  fire  allocation  Is  uniform  with  respect  to  weapon  type  and 
that  kills  within  a type  of  weapon  are  blnomlally  distributed. 
That  Is,  for  each  1 and  J, 


(l-qg(i,J))"qg(l,J)^-^  " 


for  n » 0,  ...,  yj . The  interpretation  Is  this:  when  a Blue 

weapon  of  type  1 Initiates  an  engagement,  the  type  J of  the 
weapons  to  be  engaged  is  chosen  uniformly  on  the  basis  of  the 
numbers  of  opposition  weapons  of  all  types  present;  the  number 
of  type-j  weapons  destroyed  Is  blnomlally  distributed  with 
parameters  (yj ,qg(l ,J ) ) . A physical  situation  that  may  be  com- 
patible would  Involve  ammunitions  that  are  highly  target  speci- 
fic or  target  weapons  physically  grouped  by  type.  Denote  by 
Qj^(J»i)  the  corresponding  Individual  weapon  kill  probabilities 
for  Red.  The  resulting  Markov  attrition  process  Is  an  amal- 
gamation of  Homogeneous  Process  ^ and  Heterogeneous  Process  1 
and  has  Infinitesimal  generator  A given  by 


Ill  III iKt  III  If -I-  I ^ k . M 'Itl  - aiLUl' 


r 


A((x,y) ; (x;yj^,. . . . . ,yjj^)) 


y,-n 


^ r3  r3  rd»(l)[l-(l-qg(l,J))  J] 

A((x,y);  (x,y))  • - I I x y 

i«l  J-1  ^ I y. 


+ dg(j)[l-(l-qR(J,i))  ""I 


h_- 


A((x,y);  (x^,. . . ,x^_j^,m,x^^j^,.  . . ,Xj^,^;y)) 


X.  ”3 


^ 3 /XA  „ x.-m 

— r*(J  jn  j »1) ) Qp(J»i)  yj 


These  expressions  are  valid  for  0<n<yj  and  0<m<Xj^.  Again  for 
the  sake  of  brevity,  we  omit  the  corresponding  expressions  for 
the  Jump  function  and  the  transition  matrix  of  the  embedded 
Markov  chain. 

The  family  of  processes  that  can  be  constructed  In  the  man- 
ner we  have  used  above  Is  limitless,  of  course;  the  examples 
above  serve  only  to  Illustrate  what  can  be  done  In  the  simplest 
cases.  In  Section  U we  oonslSer  the  representation  of  some 
other  specific  effects  of  Interest. 
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FURTHER  ASPECTS 


This  section  contains  additional  discussion  on  a number 
of  points  raised  In  preceding  sections. 

Computational  Problems 

In  order  to  apply  the  models  derived  here  to  computerized 
combat  simulations  one  must  develop  means  of  approximating  ex- 
pected numbers  of  survivors  at  various  times.  To  Indicate  how 
this  problem  Is  approached  within  the  context  of  the  model  of 
Section  2 without  obscuring  the  discussion  by  unnecessary 
generality,  let  us  restrict  attention  to  a homogeneous  process. 
Let  denote  the  probability  law  of  the  process  with  Ini- 
tial conditions  Bq  ■ 1,  Rq  ■ J,  the  corresponding  ex- 

pectation operator,  and  (P^)  the  transition  function  of  the 
process.  The  main  problem  Is  computation  or  approximation  of 
the  expected  number  of  Blue  weapons  surviving  at  a fixed  time 
t,  namely 

e(1*J)(;bJ  - ? k i P.((l,J),(k,0)  ; 

^ k-1  4-1  ^ 

the  approximation  of  Is  entirely  analogous.  The  given 

data  are  the  entries  of  the  Infinitesimal  generator  A.  It 
follows  from  the  backward  and  forward  equations  that  for  each  t, 

P - • I tV/n!  . (6) 

^ n-0 

This  computation  Is  usually  Impossible  to  carry  out  in  closed 
form,  but  one  may  approximate  the  Infinite  sum  by  various 
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finite  sums.  For  example,  for  each  (1,J)  and  (k,£) 

(k.O)  ~ I((l,J),(k,0)  + tA((l,j),(k,A)) 
so  that,  consequently, 

. 1 + t I k i A((l,J),(k,0)  . (7) 

^ k-0  t-o 

To  illustrate  the  meaning  of  (7),  we  list  in  Table  3 below  the 
specific  form  it  takes  for  each  of  the  four  homogeneous  pro- 
cesses presented  in  Section  3. 


Table  3.  APPROXIMATIONS  TO  E[B.] 


Process 

Expression  (7)  for  This  Process 

1 

2 

E(i.j)[B^]  - 1 . Ijd^k^t 

3 

E<^*J)[B^]  ~ 1 - Ud^k^t 

4 

E(1.j)[B^]  -V  i . Ijrgq^t 

The  approximations  generated  by  (7)  are  computationally 
useful  in  that  no  storage  requirements  and  no  long  calculations 
are  Involved.  The  same  need  not  be  true,  for  example,  of  the 
second  order  approximation 


:^i  + t I k i 


k-0  l-O 


A((l,J),(k,M) 


+ ^ i Ic  i A^((l,J),(k,D)  . 


<8) 


k-0  i-0 


28 


So  long  as  the  generator  matrix  A Is  available  In  closed  form 

(as  It  Is  In  all  the  processes  obtainable  within  the  general 

setting  of  Section  2)  the  storage  requirements  for  A Itself 

are  minimal.  But  A may  not  admit  a simple  closed  form 

expression  and  so  may  have  to  be  stored  In  a computer  In  order 

to  make  use  of  (8);  the  burden  could  be  significant.  In  addl- 

2 

tlon,  substantial  computation  time  Is  needed  to  calculate  A 
If  the  number  of  weapon  types  Is  even  moderately  large.  Finally, 
other  assumptions  may  negate  the  better  degree  of  approximation 
of  (8)  as  opposed  to  (7).  For  purposes  of  comparison  we  give 
In  Table  4 the  approximation  to  obtained  by  applying 

(8)  to  the  first  two  homogeneous  processes  described  In 
Section  3. 


Table  4.  FURTHER  APPROXIMATIONS  TO  E[B^] 


Process 

Expression  (8)  for  This  Process 

1 

2 

- 1 - J»*RPRt  + 1»'bPb'’rPr 

2 

- uw  * t 

A heuristic  Justification  of  the  first  entry  In  Table  4 
Is  the  following:  the  first  order  estimate 

e(1»J)[b^]  ^ 1 - Jrj^P^t  (9) 

in  fact  overestimates  Blue  losses,  since  Red  strength  does  not 
remain  at  J throughout  the  time  Interval  [0,t].  Using  the 
analogous  estimate 

e(1,J)[r^]  ~j  - irgPgt  , 
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one  can  take  as  a measure  of  Red  strength  the  average 


J » j[J  + (J-irgPgt)]  - J 1-5_  ; (10) 

replacement  of  J by  J In  (9)  yields  the  first  expression  in 
Table  4.  The  expression  for  Process  2 is  derived  in  an 
analogous  manner. 

The  point  is  that  from  a computational  standpoint  the  first 
order  approximation  (7)  is  always  feasible  and  is  directly  avail- 
able from  the  infinitesimal  generator  of  the  process.  Higher 
order  approximations  may  or  may  not  be  computationally  feasible; 
if  such  approximations  can  be  obtained  in  closed  form  (as  in 
Table  4)  their  use  is  probably  Justified.  Within  theater-level 
modeling  contexts  use  of  approximations  of  third  degree  or 
higher  is  almost  certainly  as  unnecessary  as  it  is  difficult. 

Another  still  largely  unexplored  computational  tool  is  the 
embedded  Markov  chain  associated  with  the  attrition  process 
under  study.  Distributions,  expectations  and  variances  of  ran- 
dom variables  defined  in  terms  of  the  attrition  process,  but  not 
explicitly  Involving  the  continuous  time  parameter,  can  be  cal-  / 
culated  using  well-developed  computational  methods  for  Markov 
chains.  For  example,  suppose  there  is  defined  a fixed  subs^  A 
of  the  state  space  E which  is  a termination  set  in  the  s^se 
that  the  engagement  terminates  at  the  random  time  T at  yhich  the 
process  enters  A.  Then  the  distribution, 

expectation,  and  variance  of  the  terminal  state  (B^,R^)  could  be 
obtained  using  Markov  chain  methods.  The  length  T of  the  engage- 
ment could  not  be  so  treated.  For  some  examples  of  such  analyses 
the  reader  Is  referred  to  [10]  and  [11]. 
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Fire  Allocation 

As  mentioned  In  Section  3»  derivation  of  the  fire  alloca- 
tion distributions  for  weapons  of  class  IS  Is  a problem 

of  considerable  practical  Interest;  In  the  heterogeneous  pro- 
cesses discussed  there  we  considered  only  the  uniform  fire 
allocation 

ag(l,yiJ)  * yj/  I y^^  . 

This  particular  form  of  allocation  Is  not  always  appropriate  so 
we  shall  briefly  discuss  some  alternatives;  these  Illustrate  the 
variety  of  effects  that  can  be  represented  within  the  context  of 
our  attrition  process.  For  the  rest  of  the  discussion  we  con- 
fine ourselves  to  a particular  type  1 of  Blue  weapons. 


1 . Restricted  Uniform  Allocation.  Suppose  there  exists  a 
subset  of  the  Index  set  of  Red  weapon  types  such 

that  a Blue  weapon  of  type  1 will  never  engage  a Red  weapon  of 
type  J If  J ^ ; such  Red  weapons  may  be  Interpreted  as  Invul- 

nerable to  Blue  weapons  of  type  1.  Otherwise,  the  fire  alloca- 
tion Is  uniform.  Then  we  would  have 


ag(l,y;J) 


y^/  I y.  If  J € K. 

J . c Tf  X.  -L 


otherwise  . 


2.  Priority  Fire  Allocation.  Suppose  that  the  N Red  weapon 
types  are  (unambiguously)  ranked  In  some  order  ^ * 

where  (J ^^(1) , . . . ,Jjj(l) ) Is  simply  a permutation  of  (1,...,N). 
Further  suppose  that  a Red  weapon  of  type  J,,(l)  will  never  be 
engaged  If  a Red  weapon  of  type  present.  This  would 

be  represented  by 

ag(l,y ; J|^n(l) ) ■ 1 , 
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where 


k*  = mln{k:  yj  (i)^0). 

k 

The  type-i  weapons  concentrate  their  fire  upon  the  highest 
priority  opposition  weapons  currently  surviving. 

3.  An  Axiomatic  Fire  Allocation.  We  note  here  that  all 
allocations  satisfying  a simple  set  of  axioms  are  of  a particu- 
lar form.  The  axioms  are; 

(1)  ag(i,y;J)  = 0 if  and  only  if  yj  = 0; 

(11)  For  all  target  forces  y and  z and  all  J, 

[N  "1 

I ap.(l,y+z;«,) 

n»i  y,  + ® J 

r N z„  1 


r N z„  1 

Uh  ^SL  * ^ J 


We  Interpret  these  axioms  in  physical  terms  as  follows.  The 
first  states  that  a weapon  type  not  present  receives  no  fire 
but  that  every  weapon  type  present  receives  a positive  fraction 
of  the  fire.  The  second  is  best  explained  step  by  step.  Con- 
sider the  effect  of  combining  two  target  forces  y and  z into 
the  single  force  w « y + z.  ag(l,y+z;Jl)  is  the  proportion  of 
fire  directed  at  the  combined  force  w which  is  allocated  to 
type  I weapons.  If  this  fire  is  further  allocated  among  the 
type  I weapons  from  tne  two  component  forces  y and  z in  propor- 
tion to  the  relative  numbers  of  such  weapons  present,  then 

— 


Is  the  fraction  of  fire  directed  at  the  combined  force  that  Is 
allocated  to  type  I weapons  originally  part  of  the  y-force. 
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Thus 

N y 

a = I y TV' 

is  the  proportion  of  fire  allocated  to  weapons  of  types  originally 
in  the  y-force  and 


a.  = I 


1=1  ^1*^1  ^ 


a^,(l,y+z;«-) 


the  fire  allocated  to  weapons  originally  in  the  z-force.  But  the 
fire  represented  by  a should  be  allocated  among  weapon  types 

y 

according  to  the  distribution  ao(l,y;*)  and  similarly  for  a so 

n z 

that  one  should  have 

ag(i,y+z;-)  = ayag(l,y;-)  + a^ag(i,z;*)  . 

If  an  assumption  of  this  form  did  not  hold,  consistency  problems 
would  arise,  with  fire  allocation  dependent  on  names  given  tar- 
gets rather  than  only  numbers  of  targets.  (For  example, 
arbitrarily  splitting  a class  of  n indistinguishable  weapons 
into  two  subclasses  of  n^  "Type  A"  and  n - n^  "Type  B"  weapons 
would  change  the  fire  allocation.) 

It  has  been  shown  by  J.  BLANKENSHIP  [1]  that  the  axioms 
(1)  and  (il)  imply  that 


Where  bj  « ag(l, (1, . . . ,1) ; J ) . In  this  case,  the  fire  allocation 
for  any  target  force  is  specified  by  that  for  any  given  base  force, 
and  is  linear  in  numbers  of  target  weapons. 
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Km  Distributions 


Analyais  of  the  sort  used  to  define  different  forms  of  fire 
allocation  Is  also  necessary  to  derive  plausible  specific  forms 
for  the  kill  distributions  for  weapons  of  classes  IM  and  PM, 
especially  the  former.  Unfortunately,  little  such  analysis  has 
been  carried  out.  In  particular,  some  dependences  allowed  by 
the  notations  can  probably  be  neglected.  For  example,  weapons 
of  class  PM  might  be  allowed  to  kill  only  weapons  of  the  same 
type  as  the  weapon  whose  detection  initiated  the  engagement. 

If  the  binomial  distributions  of  Homogeneous  Process  4 were 
felt  to  be  applicable  then  we  would  have 

for  z « (yj^. . . . • • • ,yjj)  and  4 = 0,  ...,  yj.  Further 

I work  along  these  lines  is  clearly  necessary  in  order  to  fully 

exploit  the  capabilities  of  our  model. 

Further  Phenomena 

It  is  assumed  throughout  that  all  engagements  occur  instan- 
taneously, with  ensuing  total  loss  of  contact.  A method  for 

relaxing  this  assumption,  but  retaining  the  Markov  property  of 
the  attrition  process,  is  discussed  in  [8].  In  [8]  a discussion 
is  also  given  of  possible  Inclusion  of  randomly  arriving  rein- 
forcements. Neither  of  these  problems,  however,  is  of  the  same 
Interest  or  importance  as  derivation  of  kill  distributions  for 
multiple-kill  weapons. 

A final  problem,  which  could  turn  out  to  be  of  considerable 
significance,  concerns  dependences  not  allowed  in  the  process  of 
Section  2.  None  of  the  parameters  depends  on  the  structure  of 
the  force  to  which  the  weapon  Initiating  an  engagement  belongs. 

A plausible  representation  of  such  dependence  would  permit  model- 
ing of  synergistic  relations  among  weapons  on  the  same  side,  a 
phenomenon  believed  to  be  of  considerable  importance  but  so  far 
not  well  treated  in  attrition  models. 
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